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ABSTRACT 

Aims. We investigate the ratio of the vertical and horizontal velocities of the photospheric plasma flows in the region of emerging 
magnetic fields on the Sun. 

Methods. We carried out a study of photospheric velocities during the first hours of the appearance of 83 active regions with magnetic 
flux more than 10"' Mx with data acquired by the Michelson Doppler Imager (MDI) on board the Solar and Heliospheric Observatory 
(SOHO). The emerging magnetic fluxes under investigation were isolated from extended concentrations of existing magnetic fields; 
they have different spatial scales and are located at different distances from the solar disk center. 

Results. We found that the values of maximum negative Doppler velocities that accompany the emergence of active region magnetic 
fields during the first 12 hours increase nonlinearly with the heliocentric angle. This result shows that the horizontal photospheric 
velocities of plasma outflows are higher than the vertical ones of the plasma upflows during the first hours of the emergence of active 
regions. The horizontal velocity component at the site of emerging active regions exceeds that of convective flows in the quiet Sun. A 
comparison between the velocities and the maximum value of the total magnetic flux has not revealed any relation. 
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1. Introduction 

The emerging magnetic flux of active regions on the solar sur- 
face appears in the form of separate loops; they then coalesce 
to form opposite polarity poles that separate and increase in size 
( Strous & Zwaan|[l999[ l. Earlier studies of velocities of plasma 
motions were carried out in the region of emerging magnetic 
fields at different height levels of the Sun. Investigations into the 
chromospheric level have been made by Bruzek ( 1967| [1969 1 
and others (see review by ' Chou|1993) l. Velocities at the photo- 
spheric level were measured with different methods (discussed in 
more detail below). In the last decade, helioseismological meth- 
ods have been applied to study the subphot ospheric level (|Chang| 
et al.|I999[|Kosovichev et aL|200 0, Zhark ov & Thompsori|200"8[ 



Kosovichev|2009t [Komm eTaLp008„2009F " 



Direct measurements of photospheric velocities revealed 
negative values (blue Doppler shift) on the polarity inversion 
line of emerging magnetic fields. The events under consideration 
were located in the central part of the solar disk, and, therefore, 
vertical motions were measured. Negative velocities or plasma 
upflows up to 1 km s"' were observed in separate magnetic loops 
inside emerging active regions (Tarbell et al. 1990; |Lites et al. 
T998||Sti-ous & Zwaantl999: Kubo et al.|2003,,Grigor'ev et al. 



2009 1. Grigor'ev et al.| ( 2007| i revealed high velocities of about 
1 .7 km s ' at the beginning of the powerful active region emer- 
gence NOAA 10488 at heliographic coordinates N08 E31 (Bo - 
+ 4.9). 

Indirect measurements of the horizontal velocities of the 
photospheric flow that accompanies the emergence of active 
regions were made by tracking the displacement of individual 
magnetic ele ments. T he ve locities obtained ranged from 0.1 to 



1.4 km s (Frazier 1972 Schoolman 1973 Strous & Zwaan 



1 1 999| l . [Barth & Livi| ( | 1 990| l studied 45 bipolar pairs in the emerg 



ing active region. The opposite polarity poles were separating 
from each other with velocities of 0.5-3.5 km s"', decreasing 
with time (consequently, the drift with regard to the polarity in- 
version line is only half that value). Grigor'ev et al.| ([2009 ) calcu- 
lated the separation velocities of the external boundaries of the 
photospheric magnetic flux in the active region NOAA 10488. 
The velocities decreased as the magnetic fields emerged: they 
were 2-2.5 km s ' by the end of the first hour and 0.3 km s ' in 
two hours and a half. 

It should be noted that there are several papers dealing with 



the study of vertical (|Guglielmino et al. |2006) an d horizon 



tal (Harvey & Martin 1973 Chou & Wang 1987 Hagenaar 



200 1 1) velocities at photospheric level during the appearance of 



ephemeral active regions. We do not discuss these works in de- 
tail, since they concern another spatial scale of emerging mag- 
netic fields. 

The papers above are the researches into photospheric veloc- 
ities that accompany the emergence of magnetic fields in active 
regions located only in the central zone of the solar disk. Plasma 
motions in the emerging magnetic fields distant from the disk 
center have never been studied before. 



2. Data processing and investigated objects 

We used full solar disk magnetograms and Dopplergrams in 
the photospheric line Ni I 6768 A and continuum images ob- 
tained on board the space observatory SOHO/MDI (Scherrer 



et al. II 1995 1. The temporal resolution of the magnetograms and 



Dopplergrams is 1 minute, that of the continuum, 96 minutes. 
The spatial resolution of the data is 4", the pixel size is approx- 
imately 2". Magnetograms with a 1.8.2 calibration level were 
used ( [Ulrich et al.|2009] ). Besides photospheric velocities of so- 
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Fig. 1. Active region NOAA 9456: (a) the solar surface 30 minutes before the magnetic field emergence; (b) 12 hours after the 
beginning of the emergence; (c) near the maximum of evolution. The left column represents Dopplergrams with magnetic field 
isolines +60 G (white line) and -60 G (black line); the positive velocities coiTespond to the motion of matter away from the observer, 
while the negative velocities coiTespond to the motion of matter toward the observer The right column represents magnetograms 
of the longitudinal field; the black color corresponds to the negative magnetic field, while the white one con^esponds to the positive 
magnetic field. The small rectangle marks the calculation area of the maximum negative velocity during the first 12 hours of 
the magnetic field emergence; the big one indicates the calculation area of the total magnetic flux. Plots: (d), (e) variation in the 
parameters investigated during the first 12 hours of the magnetic field emergence; (f), (g) variation in the parameters investigated 
until the maximum of the active region evolution. The beginning of the active region appearance is marked by the vertical solid line; 
the extreme values of the maximum negative velocity in plot (e) and the total magnetic flux in plot (f) are shown by the vertical 
dashed fine. 



lar plasma motions, the Dopplergrams include the contribution 
of 1) the velocity of the differential solar rotation; 2) the veloc- 
ity of the SOHO satellite with regard to the Sun; 3) the instru- 
mental distortions caused by nonuniform transmission of MDI's 
filter systems on the field of view. The technique described in 
|Grigor'ev et al.| ( [2007| l was used to separate photospheric veloc- 
ities. The negative velocity on SOHO/MDI Dopplergrams cor- 
responds to the blue Doppler shift (motion of matter toward the 
observer); the positive one, to the red Doppler shift (motion of 
matter away from the observer). There are regular errors in ve- 
locity measurements at different distances from the disk cen- 
ter For example, Doppler velocities near the limb ar e underesti- 
mated by up to 40 m s ' because of limb darkening ( Albregtsen 



& Andersen||1985| l. Additionally, we deal with different height 



levels above the solar surface in measurements of velocities at 
the disk center and near the limb because of the different optical 
thickness of the atmosphere. 

We study the velocities that accompany the emerging mag- 
netic fields. Therefore, the precise spatial superimposition of the 
used data is vary important. For this purpose, the region of the 
emerging magnetic fields was taken from the time sequences of 
magnetograms and Dopplergrams, taking into account solar ro- 
tation. An approximate value of the displacement of the region 
was calculated with the differential rotation law for photospheric 
magnetic fields (Snodg rass|1983] ). The exact tracking of the re- 
gion under investigation was performed by applying two mag- 
netograms adjacent in time, with the use of cross-correlation 
analysis. This procedure requires the existence of the magnetic 
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poles which slightly vary in time on the magnetograms. The he- 
liographic coordinates obtained were used to crop identical frag- 
ments with a size of 160"xl60" (or 80x80 pixels) from the mag- 
netogram and Dopplergram obtained at the same time. Thus, we 
achieved the precise spatial superimposition of the data. It has 
allowed us to draw reliable conclusions about the processes that 
take place. For the active region emerging near the limb, we se- 
lected the cropped region in a way that it excludes the area out- 
side the limb. 

Time variations of the total magnetic flux of the active region 
and the maximum negative velocities were calculated with the 
use of a sequence of the magnetogram and Dopplergram frag- 
ments corresponding to the emerging active region. The calcu- 
lation area was limited to the region of the emerging magnetic 
fields (Fig. [TJ. The position of the limits was visually controlled. 

The plasma motions connected with the emerging magnetic 
fields are observed on the background of the convective flows. 
If the average velocity or the velocity flux is considered as a 
characteristic of motions of matter, the center-limb dependence 
of convective velocities will make a substantial contribution. 
Besides, the active regions under consideration have different 
spatial scales and are located at different distances from the solar 
disk center; therefore the solar surface regions where the calcu- 
lation is performed have different sizes. This complicates a com- 
parison of the emerging active regions that are located at differ- 
ent distances from the disk center. The maximum negative veloc- 
ity (motion of matter toward the observer) is taken as a character- 
istic of plasma motion in the regions of the emerging magnetic 
fields during the first 12 hours (Fig. [TJd). The magnetic field 
emergence begins with the appearance of the loop apex where 
the magnetic field is horizontal; therefore, we began the calcu- 
lation of the maximum velocity values 30 minutes before the 
beginning of the appearance of active region. The magnetic flux 
of the active region emerges as separate loops, and the maximum 
negative velocity corresponds to the appearance of a single mag- 
netic loop. Therefore, the chosen parameter does not character- 
ize the emerging active region as a whole, it only shows the high- 
est velocity. For example, two peaks corresponding to the emer- 
gence of two different magnetic loops can be distinguished in the 
plot of variations in the maximum negative velocity during the 
first hours of the NOAA 9456 appearance (Fig. [T^). The residual 
magnetic flux emerges at lower velocities. To exclude Evershed 
flows we also took care of the time of the penumbra formation 
in sunspots, using continuum images. The Evershed flows are 
horizontal photospheric outflows in a sunspot penumbra. Their 
maximum velocity values can reach 2 km s ' in SOHO/MDI low 
spatial resolution data ([ Bai et al.|1998 1. 

The total magnetic flux is calculated inside isolines +60 G 
taking into account the projective effect (Fig. [T];) 



= \J](BrSi)\ + \Y,(BrSd\, 



(1) 



where O is the total magnetic flux in Mx, B, is the magnetic field 
induction in G, S , is the area on the solar surface corresponding 
to / - th pixel, «+ is the number of pixels inside isoline H-60 G 
(positive polarity) and «_ is the number of pixels inside isoline - 
60 G (negative polarity). The isoline level of 60 G was chosen to 
exclude both the magnetogram noise (about 30 G) and the con- 
tribution of the magnetic field with a short lifetime. Sometimes 
small magnetic loops appear in the place of the future active re- 
gion several hours prior to the main emergence of the magnetic 
fields. The beginning of the active region emergence was taken 
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Fig. 2. Total magnetic flux of the active regions under investiga- 
tion versus the heliocentric angle that corresponds to the begin- 
ning of the magnetic field emergence. The total magnetic flux 
maximum was determined by the inflection point in the increase 
of the flux curve (marked by diamonds) or by the last measure- 
ment (marked by asterisks). Active regions with magnetic satu- 
ration are also marked by asterisks. In active regions marked by 
asterisks, the maximum of total magnetic flux may be higher. 



as a time moment corresponding to the beginning of the con- 
tinuous growth of magnetic flux. We determined the maximum 
value of the total magnetic flux by the first inflection point of 
the curve where the flux increase was followed by its decrease 
(Fig. [T]0 or by the last value of the plot (for active regions pass- 
ing beyond the west-limb or with insufficiently downloaded data 
sequences). Note that powerful active regions significantly in- 
crease in size as they develop, and the extended concentrations 
of the magnetic fields existing on the surface enter into the re- 
gion of new flux emergence (Fig. [T^, [T]?, and [T]:). Interaction 
between existing and emerging magnetic fields is accompanied 
by integration or cancellation processes and, hence, by increase 
or decrease in the magnetic flux. In order to take this into ac- 
count to some extent, the signal background that existed before 
the emergence of magnetic fields was subtracted from the total 
magnetic flux maximum in all active regions. 

The selection of active regions for this investigation was 
based on the following criteria. Active regions have to emerge 
on the visible side of the solar disk. Emerging magnetic fields 
have to be isolated from the extended concentrations of exist- 
ing magnetic fields. The presence of single poles with a total 
magnetic flux of less than 0.5x10^' Mx was allowed in the area 
of the direct emergence of an active region during the first 12 
hours. The objects under investigation have to have complete 
data series with a one-minute time resolution during the first 12 
hours of appearance. With this in mind, the 83 active regions 
were selected for the period of 1999-2008. Among them, there 
were 80 active regions with a normal magnetic field configura- 
tion and three with an inverse magnetic field configuration ac- 
cording to Hale's polarity law. The events had different spatial 
scales and were located at different distances from the solar disk 
center (Fig. |2]i. The total magnetic flux of each event exceeded 
10^' Mx. Magnetic saturation is possible in SOHO/MDI mea- 
surements of the magnetic field strength. This occurs when the 
spectral line is shifted out of the filter system passband. |Liu et alT] 
(2007) modeled MDI measurements and demonstrated that mag- 
netic saturation takes place only in strong magnetic fields with 
high velocities. According to their calculations, magnetic satu- 
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Fig. 3. (a) Maximum negative velocities during the first 12 hours of the emergence of active regions versus the heliocentric angle 
that corresponds to the position of the active regions at this point in time. Diamonds mark normal active regions; black squares mark 
inverse ones. The solid line corresponds to the nonlinear regress equation (2); dotted lines correspond to 99% confidence intervals 
for the means, (b) Maximum negative velocities of convective flows in the quiet Sun versus the heliocentric angle. 



ration is reached at the magnetic field strength of 2800 G, when 
the velocity is +2000 m s"'. In our study, the maximum mag- 
netic field strength exceeds 2800 G only in two events of the ac- 
tive regions under consideration. These events are also marked 
by asterisks in Fig. [2] The morphology of the active regions 
(i.e., the degree of fragmentation of magnetic fields, the rate and 
monotony of magnetic flux increase, and the orientation of the 
axis connecting opposite polarity poles with regard to the line of 
sight) was not taken into account. 



3. Results 

Figure [3^ presents the maximum negative velocities obseiTed 
during the first 12 hours emergence of active regions according 
to their position with regard to the disk center. The distance from 
the disk center is expressed in the heliocentric angle 6* - it is the 
angle between the normal to the surface and the line-of-sight 
to the emerging magnetic flux. The maximum negative veloci- 
ties obtained are observed, as a rule, on the polarity inversion 
line of the emerging magnetic fields; however, they are related 
to the convection in adjacent regions for some events in the cen- 
tral part of the disk (0 < 25°). Figure [3^ shows that values of 
maximum negative velocities increase nonlinearly with the he- 
liocentric angle. The mean velocity for - 20° is -673 m s"'; 
the deviation from the mean does not exceed 200 m s The 
mean velocity for - 60° is -1277 m s"'; the deviation from 
the mean reaches up to 550 m s One can see a growth ten- 
dency of deviation from mean with increasing (Fig. [3^). The 
high deviation for large indicates the existence of other rela- 
tions that determe the value of the horizontal velocity of plasma 
flow at the site of emerging magnetic fields at the photospheric 
level. Our statistics include three inverse active regions whose 
parameters do not deviate from the dependence formed by nor- 
mal active regions. Thus we see that the horizontal velocities of 
plasma outflows exceed the vertical ones of plasma upflows at 
the beginning of the emergence of active regions. 

Figure ^ shows the maximum negative velocities for the 
quiet Sun connected with convective flows. To calculate this de- 
pendence, the region on the surface was traced taking into ac- 



count solar rotation. The region with a size of 40" X 40" under 
study was going through the longitude range of W00-W67 at the 
latitude of S04 from 3 to 8 May 1999. The tilt of the solar North 
rotational axis toward the observer Bq was -3.5°. Data with a 
one-minute resolution were used; they were processed with the 
technique described in section 2 of this paper. In Fig. ^ one can 
see that the velocities increase with the heliocentric angle; these 
values do not exceed 1200 m s"'. The comparison between the 
plots in Fig. [3^ and Fig. ^ has revealed that the horizontal ve- 
locity component at the sites of the emergence of active regions 
is higher than that of the convective flow in the quiet Sun. 

We performed the regression analysis with the second degree 
polynom. As a result, we obtained the equation 

V = 6.01 -40.180 + 0.3 16l^ (2) 

where v is the maximum negative Doppler velocity during the 
first 12 hours of the emergence of active regions in ms"' and 
is the heliocentric angle in degrees. According to the F-statistic 
this regression model is significant. The correlation ratio (es- 
timate of the closeness of the nonlinear relation) between the 
value of the maximum negative velocity and the heliocentric an- 
gle 6 (or the angle of view to the emerging magnetic flux) is 
equal to -0.74; this implies a high relation between the parame- 
ters under consideration. The confidence intervals for the means 
were calculated with a confidence probability of 99%. They are 
marked by dotted lines in Fig. [3^. The confidence intervals have 
close limits for 10° < 9 < 70° with high statistics and show that 
equation (2) approximates the data well. The limits of the con- 
fidence intervals essentially increase for < 10° and 70° < 
with low statistics. It follows from equation (2) that the mean 
vertical velocity for = 0° is 6 m s"' with a confidence interval 
of ±399 ms"^ Positive velocity values are false in this confi- 
dence range. Obviously, velocities coiTesponding to = 0° at 
the sites of emerging active regions will be not lower than the 
velocities of convective flows of the quiet Sun ( — 150 m s"'). 
The absence of statistics for < 10° at the emerging active re- 
gions does not allow one to take it into account in regression 
equation (2). The mean horizontal velocity for = 90° from 
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Fig. 4. Maximum negative velocities during the first 12 hours of 
the magnetic field emergence versus the maximum value of the 
total magnetic flux in the active regions. The designations are the 
same as for the Fig. |2j 



equation (2) is -1081.44 m s"' with a confidence interval ±601 
m s ' . The wide confidential interval shows that this value of the 
mean horizontal velocity is ambiguous. 

A comparison of the maximum negative velocities with the 
total magnetic flux did not reveal any relation (Fig. |4|i. The com- 
parison between the plots in Fig. |2]and Fig. [3^ also showed no 
connection, because high velocities were not observed in emerg- 
ing active regions with high magnetic flux in the central part of 
the solar disk, whereas near the limb there may be high velocity 
values even in active regions with low magnetic flux. From the 
author's point of view the absence of a dependence between ve- 
locity and total magnetic flux can be explained by the fact that 
active region magnetic fields emerge by separate fragments on 
time scales from several hours to 5-7 days and therefore the 
maximum negative velocity during the first 12 hours does not 
characterize an active region emergence as a whole. 



4. Conclusions 

We presented a statistical study of the velocities of plasma flows 
that accompany the emergence of active regions during the first 
12 hours at the photospheric level with data with high tempo- 
ral resolution. We found that the values of maximum negative 
Doppler velocities increase nonlinearly with the heliocentric an- 
gle. This shows that the horizontal velocities of plasma outflows 
exceed the vertical ones of plasma upflows during the first hours 
of active region emergence. Horizontal velocities at the sites of 
emerging active regions are higher than those of convective flows 
in the quiet Sun. This result is a direct confirmation of theoreti- 
cal models that showed that magnetic fields emerging in the solar 
atmosphere expand faster in the horizontal direction than in the 
vertical direction (e.g. [Shibata et al.|l989[|Fan|200l||Magara & 
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